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ABSTRACT 

We present a measurement of the cosmic microwave background (CMB) temperature power spec- 
trum using data from the recently completed South Pole Telescope Sunyaev-Zel'dovich (SPT-SZ) 
survey. This measurement is made from observations of 2540 deg 2 of sky with arcminute resolution 
at 150 GHz, and improves upon previous measurements using the SPT by tripling the sky area. We 
report CMB temperature anisotropy power over the multipole range 650 < I < 3000. We fit the 
SPT bandpowers, combined with the results from the seven-year Wilkinson Microwave Anisotropy 
Probe ( WMAP7) data release, with a six-parameter ACDM cosmological model and find that the two 
datasets are consistent and well fit by the model. Adding SPT measurements significantly improves 
ACDM parameter constraints, and in particular tightens the constraint on the angular sound horizon 
9 S by a factor of 2.7. The impact of gravitational lensing on the CMB power spectrum is detected with 
8.1 fx, the most significant detection to date. The inferred amplitude of the lensing spectrum is consis- 
tent with the ACDM prediction. This sensitivity of the SPT+ WMAP7 data to lensing by large-scale 
structure at low redshifts allows us to constrain the mean curvature of the observable universe with 
CMB data alone to be Q k = -0.003^0^. Using the SPT+WMAP7 data, we measure the spectral 
index of scalar fluctuations to be n s = 0.9623 ± 0.0097 in the ACDM model, a 3.9 a preference for a 
scale-dependent spectrum with n s < 1. The SPT measurement of the CMB damping tail helps break 
the degeneracy that exists between the tensor-to-scalar ratio r and n s in large-scale CMB measure- 
ments, leading to an upper limit of r < 0.18 (95% C.L.) in the ACDM+r model. Adding low-redshift 
measurements of the Hubble constant (Hq) and the baryon acoustic oscillation (BAO) feature to the 
SPT+ WMAP7 data leads to further improvements. The combination of SPT+ WMAP7+H +BAO 
constrains n s = 0.9538 ± 0.0081 in the ACDM model, a 5.7 a detection of n s < 1, and places an 
upper limit of r < 0.11 (95% C.L.) in the ACDM+r model. These new constraints on n s and r have 
significant implications for our understanding of inflation, which we discuss in the context of selected 
single-field inflation models. 

Subject headings: cosmology - cosmology:cosmic microwave background - cosmology: observations - 
large-scale structure of universe 
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1. INTRODUCTION 

Over the past two decades, measurements of the cos- 
mic microwave background (CMB) have provided pro- 
found insight into the nature of the universe. Detailed 
information about the composition and evolution of the 
universe is encoded in the temperature and polariza- 
tion anisotropy of the CMB. Measuring this anisotropy 
enables powerful tests of cosmological theory. On de- 
gree scales, CMB anisotropy is primarily generated by 
the acoustic oscillations of the primordial plasma in the 
early universe. The Wilkinson Microwave Anisotropy 
Probe ( WMAP) satellite has been used to measure these 
acoustic oscillations with cosmic variance-limited p reci- 
sion on angular s cales corresponding to I < 500 (Ko- 
matsu et al.|2011 hereafter WMAP7). On much smaller 
angular scales, primary CMB anisotropy becomes dom- 
inated by effects imprinted on the CMB at low red- 
shift (so-called secondary anisotropy) and foregrounds; 
at millimeter wavelengths, this transition occurs at £ ~ 
3000. This small-angular scale millimeter-wavelength 
anisotropy h as been measured by the South Pole Tele- 



scope (SPT) jLueker et al. 2010 Shirokoff et al. 



ichardt et ai.[ | 2012a[) and the Atacama Cosmology Tele- 



2011 



Re- 



Das et al. 



2011a 



scope (ACT) flFowler et al.l|2010| 

The anisotropy in the ClVlB at intermediate angular 
scales, 1000 < £ < 3000, is often referred to as the 
"damping tail" since the anisotropy power on these an- 
gular scal es is damp ed by photon diffusion during recom- 
bination ( Silk 1968 ) . Adding measurements of the damp- 
ing tail to large-scale CMB measurements considerably 
tightens the resulting cosmological constraints. There- 
fore measurement on smaller angular scales are sensitive 
to the photon diffusion scale during recombination. The 
wider range of angular scales also enables better con- 
straints on the sound horizon at recombination (by mea- 
suring more acoustic peaks) and the slope of the primor- 
dial power spectrum. Finally, although tensor pertur- 
bations from cosmic inflation add CMB power only at 
very large angular scales, the effect of these tensor per- 
turbations is degenerate with changes in n s in large-scale 
measurements. Damping tail measurements help break 
this degeneracy, thus tightening constraints on the level 
of tensor perturbations. 

In the past few years, there have been several in- 
creasingly precise measu rements of the CMB damping 
tail, including ACBAR (iReichardt et al. || 2009[ ), Q UaD 
dBrown et aLl|2009[ |Friedman et al.||200ifc iroT ((Das 
et al.||2011b[ ), and SPT QKeisler et al.||2011[ hereafter 
Kll). The most precise published measurement of the 
CMB damping tail prior to this work comes from the first 
790 deg 2 of the South Pole Telescope Sunyaev-Zel'dovich 
(SPT-SZ) survey (Kll). 

In this paper, we present a measurement of the power 
spectrum from the third acoustic peak through the CMB 
damping tail, covering the range of angular scales corre- 
sponding to multipoles 650 < £ < 3000. This power 
spectrum is calculated from the complete SPT-SZ sur- 

22 Harvard-Smithsonian Center for Astrophysics, 60 Garden 
Street, Cambridge, MA, USA 02138 

23 Berkeley Center for Cosmological Physics, Department of 
Physics, University of California, and Lawrence Berkeley National 
Laboratory, Berkeley, CA, USA 94720 



vey covering 2540 deg of sky, and improves upon the 
results presented in Kll by expanding the sky coverage 
by a factor of three. 

We present constraints from this measurement on the 
standard ACDM model of cosmology, then extend the 
model to quantify the amplitude of gravitational lens- 
ing of the CMB. We use this sensitivity to gravitational 
lensing by large-scale structure at low redshifts to mea- 
sure the mean curvature of the observable universe from 
CMB data alone. We also consider models including ten- 
sor perturbations, and explore implications of the result- 
ing parameter constraints for simple models of inflation. 
Adding low-redshift measurements of the Hubble con- 
stant (Ho) and the baryon acoustic oscillation (BAO) 
feature to the CMB data further tightens parameter con- 
straints, and we present combined parameter constraints 
for each of the above model extensions. The implications 
of the SPT power spectrum for a larger range of exten- 
sions to the standard cosmologi cal model are exp lored 
more fully in a companion paper, Hou et al. (2012 here- 
after H12). 

This paper is organized as follows. We describe the 
SPT observations and data reduction in Sj2] We present 
the power spectrum calculation in <|3| We discuss tests 
for systematic errors injp] We present the power spec- 
trum measurement in |5] In fj6j we outline our cosmo- 
logical parameter fitting framework and present the re- 
sulting parameter constraints, then use these constraints 
to explore the implications for simple models of inflation. 
Finally, we conclude in <|7j 

2. OBSERVATIONS AND DATA REDUCTION 

The SPT is a 10-meter diameter telescope located at 
the Amundsen-Scott South Pole station in Antarctica. 
The first survey with the SPT, referred to as the "SPT- 
SZ" survey, was completed in November 2011 and cov- 
ered a ~ 2500 deg 2 region of sky between declinations of 
-40° and -65° and right ascensions of 20h and 7h. The 
SPT-SZ survey is shown in Figure[l] Here we present the 
first power spectrum measurement that uses data from 
the complete SPT-SZ survey. We use data from 2540 
deg 2 of sky in this analysis. 

This work uses observations and data reduction meth- 
ods that are very similar to those described in Kll. In 
this section, we give an overview of the observations 
and data reduction, highlighting the differences with the 
treatment in Kll, to which we refer the reader for a de- 
tailed treatment of the analysis methods. 

2.1. Observing Strategy and Fields 

From 2008 - 2011, the SPT was used to observe a con- 
tiguous ~ 2500 deg 2 patch of sky to a noise level of ap- 
proximately 18 /iK-arcmin 24 at 150 GHz. 25 This area 
of sky was observed in 19 contiguous sub-regions which 

24 Throughout this work, the unit K refers to equivalent fluctu- 
ations in the CMB temperature, i.e., the temperature fluctuation 
of a 2.73 K blackbody that would be required to produce the same 
power fluctuation. The conversion factor is given by the derivative 
of the blackbody spectrum evaluated at 2.73 K. 

25 The SPT-SZ survey also includes data at 95 and 220 GHz. 
However, this work uses only 150 GHz data since this observing 
band is the most sensitive for the SPT and the data from one 
observing band are sufficient to make high signal-to-noise maps of 
the CMB anisotropy. 



Fig. 1. — The 2500 deg 2 SPT-SZ survey. We show the full survey region with lightly filtered 95 GHz data from the SPT, using the data 
and filters which best capture the degree-scale anisotropy of the CMB visible in this figure. The power spectrum measureme nt r eported in 
this paper is calculated from 2540 deg 2 of sky and analyzes 150 GHz data with a different high-pass filter, as described in § |2,2[ 



we refer to as observation "fields". In the basic survey 
strategy, the SPT was used to observe a single field until 
the desired noise level was reached before moving on to 
the next field. Two fields were observed in 2008, three in 
2009, five in 2010, and nine in 2011. All nine fields from 
2011 were observed to partial depth in 2010 in order to 
search for massive galaxy clusters, then re-observed in 
2011 to achieve nominal noise levels. The results of that 
bright cluster search were published in [Williamso n et alT] 
(2011). In terms of sky area, this equates to observing 



167 deg 2 in 2008, 574 deg 2 in 2009, 732 deg 2 in 2010, and 
1067 deg 2 in 2011. The fields are shown in Figure |2[ a nd 
the field locations and sizes are presented in Table [1] 

Both fields from 2008 (ra5h30dec-55 and 
ra23h30dec-55) were re-observed in later years 
to achieve lower than normal noise levels. In this 
analysis, we use data from only one year for each field 
because the beam and noise properties vary slightly 
between years. This choice simplifies the analysis with- 
out affecting the results as the bandpow er u ncertainties 
remain sample variance dominated (see ^3.5). 



The SPT is used to observe each field in the following 
manner. The telescope starts in one corner of the ob- 
servation field, slews back and forth across the azimuth 
range of the field, and then executes a step in elevation, 
repeating this pattern until the entire field has been cov- 
ered. This constitutes a single observation of the field, 
and takes from 30 minutes to a few hours, depending on 
the specific field being observed. Azimuthal scan speeds 
vary between fields, ranging from 0.25 to 0.42 degrees per 
second on the sky. The starting elevation positions of the 
telescope are dithered by between 0.3' and 1.08' to en- 
sure uniform coverage of the region in the final coadded 
map. 



In four of the 2008 and 2009 fields, ra23h30dec-55, 
RA21HDEC-60, ra3h30dec-60, and RA21HDEC-50, ob- 
servations were conducted with a "lead-trail" strategy. 
In this observation strategy, the field is divided into two 
halves in right ascension. The "lead" half is observed 
first, immediately followed by the "trail" half in a man- 
ner such that both halves are observed over the same 
azimuthal range. If necessary, the lead-trail data could 
be analyzed in a way that cancels ground pickup. In this 
analysis, we combine lead-trail pairs into single maps, 
and verify that contamination from ground pickup is neg- 
ligible - see below and §[4] for details. 

We apply several (often redundant) data quality cuts 
on individual observations using the following criteria: 
map noise, noise-based bolometer weight, the product of 
median bolometer weight with map noise, and the sum of 
bolometer weights over the full map. For these cuts, we 
remove outliers both above and below the median value 
for each field. We do not use observations that are flagged 
by one or more of these cuts. We also flag observations 
with only partial field coverage. Finally, we cut maps 
that were made from observations in azimuth ranges that 
could be more susceptible to ground pickup over the an- 
gular scales of interest. We use "ground-centered" maps 
to measure ground pickup on large [l ~ 50) scales, and 
cut observations that were made at the azimuths with the 
worst 5% ground pickup to minimize potential ground- 
pickup on smaller angular scales. Although this cut does 
have an impact on our null tests (see § Pp) , we emphasize 
that it does not significantly change the power spectrum, 
the precision of which is limited by sample variance. 

2.2. Map-making: Time Ordered Data to Maps 

As the SPT scans across the sky, the response of each 
detector is recorded as time-ordered data (TOD). These 
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TABLE 1 

The fields observed with the SPT between 2008 and 2011 



Name 


R.A. O 


Decl. (°) 


AR.A. (°) 


ADecl. (°) 


Effective Area (deg 2 ) 


ra5h30dec-55 


82.7 


-55.0 


15 


10 


84 


ra23h30dec-55 


352.5 


-55.0 


15 


10 


83 


ra21hdec-60 


315.0 


-60.0 


30 


10 


155 


ra3h30dec-60 


52.5 


-60.0 


45 


10 


227 


ra21hdec-50 


315.0 


-50.0 


30 


10 


192 


ra4h10dec-50 


62.5 


-50.0 


25 


10 


156 


ra0h50dec-50 


12.5 


-50.0 


25 


10 


157 


ra2ii30dec-50 


37.5 


-50.0 


25 


10 


157 


ra1hdec-60 


15.0 


-60.0 


30 


10 


152 


ra5i-i30dec-45 


82.5 


-45.0 


15 


10 


109 


ra6ii30dec-55 


97.5 


-55.0 


15 


10 


85 


RA23HDEC-62.5 


345.0 


-62.5 


30 


5 


75 


RA21HDEC-42.5 


315.0 


-42.5 


30 


5 


121 


ra22h30dec-55 


337.5 


-55.0 


15 


10 


84 


ra23hdec-45 


345.0 


-45.0 


30 


10 


217 


ra6hdec-62.5 


90.0 


-62.5 


30 


5 


75 


ra3h30dec-42.5 


52.5 


-42.5 


45 


5 


179 


RAlHDEC-42.5 


15.0 


-42.5 


30 


5 


119 


ra6h30dec-45 


97.5 


-45.0 


15 


10 


111 


Total 










2540 



Note. — The locations and sizes of the fields observed by the SPT between 2008 and 2011. For each field we give the center 
of the field in Right Ascension (R.A.) and declination (decl .), th e nominal extent of the field in Right Ascension and declination, 



and the effective field area as defined by the window (see £ 3 . 2 1 




Fig. 2.— The SPT was used to observe 2500 deg 2 over 19 indi- 
vidual fields, which are overlaid here on an orthographic pr ojection 
of the IRAS 100 /im dust map from |Schlegel et al.J fl998| ). These 
observation fields were chosen to lie in regions ot low dust emission. 

TOD are recorded at 100 Hz and have a Nyquist fre- 
quency of 50 Hz, which corresponds to a multipole num- 
ber parallel to the scan direction (£ x ) between 72,000 and 



43,000 at the SPT scan speeds. Since we only report the 
power spectrum up to £ = 3000, we can benefit compu- 
tationally by reducing the sampling rate. We choose a 
low-pass filter and down-sampling factor based on each 
field's scan speed such that they affect approximately 
the same angular scales. We use a down-sampling factor 
of 6 for 2008 and 2009, and 4 for 2010 and 2011 with 
associated low-pass filter frequencies of 7.5 and 11.4 Hz 
respectively. These filtering choices remove a negligible 
amount of power in the signal band. 

Next, the down-sampled TOD are bandpass filtered be- 
tween i x = 270 and 6600. The low-pass filter is necessary 
to avoid aliasing high-frequency noise to lower frequen- 
cies during map-making. The high-pass filter reduces 
low-frequency noise from the atmosphere and instrumen- 
tal readout. The high-pass filter is implemented by fit- 
ting each bolometer's TOD (from a single azimuthal scan 
across the field) to a model consisting of low-frequency 
sines and cosines and a fifth-order polynomial. The best- 
fit model is then subtracted from the TOD. During the 
filtering, we mask regions of sky within 5 arcminutes of 
point sources with fluxes of Si5o GHZ > 50 mJy. These 
regio ns a re also masked in the power spectrum analysis, 
see § |3.2| 

At this stage, the TOD retain signal from the at- 
mosphere that is correlated between detectors. We re- 
move the correlated signal by subtracting the mean sig- 
nal across each detector module for every time sample. 26 
This process acts as an approximately isotropic high-pass 
filter. 

The filtered TOD are made into maps using the pro- 

26 The SPT-SZ focal plane has a hexagonal geometry with six 
triangular bolometer modules, each with ~160 detectors. Each 
module is configured with a set of filters that determines its ob- 
serving frequency of 95, 150, or 220 GHz. 
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cess described by Kll. The data from each detector re- 
ceive a weight based on the power spectral density of 
that detector's calibrated TOD in the 1-3 Hz band. This 
band corresponds approximately to the signal band of 
this analysis. The detector data are binned into maps 
with 1' pixels based on the telescope pointing informa- 
tion. In the power spectrum analysis presented in <|3j we 
adopt the flat-sky approximation, where the wavcnum- 
ber k is equivalent to multipole moment £ and spherical 
harmonic transforms are replaced by Fourier transforms. 
We project from the curved celestial sky to flat-sky maps 
with the oblique L ambert equal-area azimuthal projec- 
tion ( |Snyder||1987| ). 



2.3. Beam Functions 

A precise measurement of the SPT beam — the opti- 
cal response as a function of angle — is needed to cali- 
brate the angular power spectrum as a function of mul- 
tipole. We summarize the method used t o measure the 
SPT b eams and refer the reader to Kll or lSchaffer et al.l 
(2011) for a more detailed description. 

The average 150 GHz beam is measured for each year 
using a combination of maps from Jupiter, Venus, and 
the 18 brightest point sources in the CMB fields. The 
maps of Jupiter are used to measure the beam outside a 
radius of 4', while the maps of the bright point sources 
are used to measure the beam inside that radius. Maps 
of Venus are used to join the inner and outer beam maps 
into a composite beam map. The maps of the planets 
are not used to estimate the very inner beam due to a 
non-linear detector response when directly viewing bright 
sources and the non-negligible angular size of the plan- 
ets compared to the beam. We use the composite beam 
map to measure the beam function Be, defined as the az- 
imuthally averaged Fourier transform of the beam map. 
We consider uncertainties in the measurement of Bi aris- 
ing from several statistical and systematic effects, such 
as residual atmospheric noise in the maps of Venus and 
Jupiter, and account for known inter-year correlations 
of some of these sources of uncertainty. The parameter 
constraints quoted in this work are not sensitive to the 
calculated beam uncertainties; we have tested increas- 
ing our beam uncertainties by a factor of two and have 
seen no significant impact on the resulting cosmological 
parameter fits. 

A nearly identical beam treatment was used by Kll. 
The main difference is that the beam function was nor- 
malized to unity at £ = 350 in Kll rather than £ = 750 
in this work. The average multipole of our calibration 
region is close to £ = 750, and this choice of normal- 
ization scale better decouples the beam and calibration 
uncertainties. 

2.4. Calibration 

The observation-to-observation relative calibration of 
the TOD is determined from repeated measu rements of 
a galactic H I I region, RCW38. As in Kll and Reichardt 



et al. (2012a 



, the absolute calibration is determined by 
comparing the SPT and WMAP7 power in several £- 
bins over the multipole range £ € [650,1000]. We use 
the same ^-bins for both experiments: seven bins with 
S£ — 50. This calibration method is model- independent, 
requiring only that the CMB power in the SPT fields 



is statistically representative of the all-sky power. We 
estimate the uncertainty in the SPT power calibration 
to be 2.6%. We do not find the parameter constraints 
quoted in this work to be sensitive to the calibration 
uncertainty; changing the calibration uncertainty by a 
factor of two in either direction has no significant impact 
on the resulting cosmological parameter fits. 

3. POWER SPECTRUM 

In this section, we describe the power spectrum calcu- 
lation. This analysis closely follows the analysis devel- 
oped by Lueker et al. (2010) and used by Kll; we refer 
the reader to those papers for a more detailed descrip- 
tion. We refer to the average power over a given range 
of £ values as a bandpower. We use a pseudo-CV method 



under the flat-sky approximation as described in j 2.2 
The power spectra are calculated independently in each 
of the 19 fields, then combined into the final result. We 
report bandpowers in terms of Dp, which is defined as 



De 



£(£ + !) 
2?r 



Ct 



(1) 



To calculate Dg, we use a cros s-spectrum bandpower es- 
timator as described in § |3.3[ which has th e adva ntage 
of being free of noise bias; see Lueker et al. (2010) for a 
more detailed description. 

3.1. Maps 

The basic input to the cross-spectrum estimator is a 
set of maps for a given field, each with independent noise. 
For most fields, this input set is comprised of maps from 
single observations. Each observation has statistically 
independent noise because observations are temporally 
separated by at least an hour and the TOD have been 
high-pass filtered at ~ 0.2 Hz. As in Kll, for the four 
fields observed with a lead-trail strategy, we construct 
the input map set by combining lead-trail pairs into sin- 
gle maps. The ra23h30dec-55 field was observed using 
comparatively large elevation steps and hence has less 
uniform coverage. For this field, the single maps that 
are the basic input to the cross-spectrum estimator are 
formed by combining two pairs of lead-trail observations. 
Each pair is chosen to have different elevation dithers, 
leading to a more homogeneous field coverage. 

3.2. Window 

For a given field, each of the maps is multiplied by the 
same window W in order to avoid sharp edges at map 
boundaries, control overlap between adjacent fields, and 
remove bright point sources. Each window is the prod- 
uct of an apodization mask with a point source mask. 
The apodization masks are calculated by applying a 1° 
taper using a Hann function to the edges of the uniform 
coverage region of each field. The observations were de- 
signed such that the uniform coverage region overlaps 
between neighboring fields. We define our apodization 
windows such that the overlap region between adjacent 
fields contains a combined weight that approaches but 
never exceeds unity (the weight at the center of the field) . 
This process results in apodization windows that include 
marginally smaller regions of sky and fall off more slowly 
than the windows used in Kll, which did not need to 
account for field overlap. 
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As was done in Kll, we identify all point sources with 
150 GHz flux > 50 mJy. Each of these point sources 
is masked with a 5'-radius disk that is tapered out- 
side the disk using a Gaussian taper with a width of 
CT taper = ^' ■ Point source masks remove 1.4% of the total 
sky area. Using previous measurements of the mm- wave 
point sourc e population (Vieira et al. 2010 Shirokoff 
et al. 2011), we estimate that the power from residual 
point sources below this flux cut is Ci ~ 1.3 X 10~ 5 //K 2 , 

or De 18 /iK 2 (5^5) ■ This power is approximately 
half the CMB anisotropy power at £ = 3000, the up- 
per edge of the multipole range reported in this analysis. 
Further discussion of the point source model is reserved 
for 6 IO 



3.3. Cross-Spectra 

The next step in calculating the power spectrum is to 
cross-correlate single maps from different observations of 
the same field. Each map is multiplied by the window 
for its field, zero-padded to the same size for all fields, 
then the Fourier transform of the map m A is calculated, 
where A is the observation index. The resulting Fourier- 
space maps have pixels of size df = 5 on a side. We 
calculate the average cross spectrum between the maps 
of two observations A and B within an ^-bin b: 



D 



AB _ 



1(1+1) 

2tt 



HftRe[rh t 



A~B* 



(2) 



■eb 



where Hi is a two-dimensional weight array described 
below, and £ is a vector in two-dimensional £-space. Each 
field typi cally has about 200 single maps in the input 
set (see § 3.1), resulting in ~ 20,000 cross-spectra. We 



average all cross-spectra D AB for A ^ B to calculate a 

binned power spectrum Db for each field. 

Due to our observation strategy, the maps have statis- 
tically anisotropic noise; at fixed £, modes that oscillate 
perpendicular to the scan direction (£ x — 0) are noisier 
than modes that oscillate parallel to the scan direction. 
This anisotropic nois e - a nd the filtering we apply to re- 
duce the noise (see § 2.2) - causes different modes in a 
given i bin to have different noise properties. As in Kll, 
we use a two-dimensional weight Hi which accounts for 
the anisotropic noise in the maps. We define the weight 
array according to 



fl/oc {CT + N t ) 



(3) 



where Cf 1 is the theoretical power spectrum used in 



simulations described in £ 3.4.1 and Ng is the two 



dimensional calibrated, beam-deconvolved noise power, 
which is calculated from difference maps in which the 
right-going scans are subtracted from the left-going 
scans. The weight array is then smoothed with a Gaus- 
sian kernel of width ag = 450 to reduce the scatter in the 
noise power estimate, and normalized to the maximum 
value in each annulus. Hi is calculated independently for 
each observation field. 

3.4. Unbiased Spectra 

The power Db is a biased estimate of the true sky 
power, Dt, due to effects such as TOD filtering, pro- 
jection effects, and mode mixing from the window. The 



biased and unbiased estimates are related by 
Db> = Kb'bDb , 



(4) 



where the K matrix accounts for the effects of the 
beams, TOD filtering, pixelization, windowing, and 
band-averaging. K can be expanded as 



K w =Pbt (M u ,[W]F e ,Bl)Q 



t'V 



(5) 



Qf/h' is the binning o perator and Pbi is its recipro- 



cal (Hivon et al 
Mi 



2002). The "mode-coupling matrix" 
[W] accounts tor mixing modes between multipole 



moments which arises from observing a finite portion of 
the sky. We calculate Mgg> [W] analytically from the win- 
dow function W foll owing the prescription described by 
Lueker et al. (2010). Over the range of multipoles re- 
ported in this analysis, the elements of the mode-coupling 
matrix depend only on the distance from the diagonal. 
Fi is the transfer function due to TOD fi ltering and map 
pixelization, which is described in { 



function described in 



3.4.1 



B\ is the beam 
For sufficiently large ^-bins, 
the K matrix is invertible, allowing an unbiased estimate 
of the true sky power: 



D h = K 



1 bb- 



,D b 



(6) 



3.4.1. Simulations and the Transfer Function 



The transfer function Fg is calculated from end-to- 
end simulations. One hundred full-sky realizations 
are generated at a Healpix 27 resolution of N s jd c =8192. 
These simulated skies include gravitationally lensed 
CMB anisotropy based on the best-fit ACDM model from 
Kll, a Poisson distribution of radio galaxies, and Gaus- 
sian realizations of the thermal and kinetic Sunyaev- 
Zel'dovich (SZ) effects and cosmic infrared background 
(CIB). The lensed realizations of the CMB spectrum ar e 
generated out to £ = 8000 using LensPix ( |Lewis| [2005 ) . 
The Poisson radio galaxy contribution is based on the |De| 
Zotti et al. ( 2010 ) model for sources below the 5 a detec- 
tion th reshold in the SPT -SZ survey, and the observed 
counts ( Vieira et al.|2010 1 above that flux. The sh ape of 
the thermal SZ spectrum is tak en from Shaw et al. ( 2010 1 
with an amplitude taken from Reichardt et al. "j 2012ap . 
T he kinetic SZ spect rum is based on the fiducial model 



Zahn et al. (2012) . The CIB spectrum is taken from 



the best-lit values m Reichardt et al. (2012a) 

Unlike the simulations in Kll, these simulations cover 
the full sky. The full-sky simulations make it simple to 
account for overlap between fields when calculating the 
sample va rian ce term of the bandpower covariance ma- 
trix (see £3.5). These simulations also account for any 



effects due to projecting from the curved sky to flat sky 
maps to first order in the transfer function, though these 
effects should be negligible, as argued in Kll. 

These simulated skies are observed using the SPT 
pointing information and then filtered and processed into 
maps using the same pipeline as for the real data. For 
each field, we calculate the transfer function by compar- 
ing the average power spectrum of these simulated maps 
to the known inpu t spectrum using an iterative scheme 
( |Hivon et aL|2002| ). 



http:/ /hcalpix.jpl. nasa.gov 
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The transfer function is equal to ~ 0.25 at £ — 650 
and reaches a plateau for I > 1200. The transfer func- 
tion does not reach unity at any scale due to the strong 
filtering of £ x < 300 modes. 

3.5. Bandpower Covariance Matrix 

The bandpower covariance matrix quantifies the bin- 
to-bin covariance of the unbiased spectrum. The covari- 
ance matrix contains signal and noise terms as well as 
terms accounting for beam and calibration uncertainties. 
The signal term, often referred to as "sample variance" , 
is c alculat ed from the 100 simulations described in sec- 
tion 3.4.1 For each simulated Healpix sky, we calculate 
the combined power spectrum from all fields, then mea- 
sure the variance of these 100 estimates. This process 
naturally accounts for any overlap between fields. The 
noise term, or "noise variance" , is estimated directly from 
the data using the dist ribution of ind i vidua l cross-spectra 



D<- as described by Lueker et al. (20101. The sample 



variance is dominant at multipoles below i < 2900. At 
smaller angular scales, the noise variance dominates. 

The initial estimate of the bandpower covariance ma- 
trix has low signal-to-noise on the off-diagonal elements. 
As in Kll, we condition the covariance matrix to reduce 
the impact of this uncertainty. 

We must also account for the bin-to-bin covariance due 
to the uncertainties in the beam function Bg. We con- 
struct a "beam correlation matrix" for each source of 
beam uncertainty: 



beam 
Pij 



so, 

D 4 



where 



(7) 



(8) 



We sum these matrices to find the full beam correlation 
matrix, and convert to a covariance matrix according to 



r-^beam beam 

*•'<.; ~ Pij 



DiDj. 



(9) 



3.6. Combining Fields 

The analysis described in the previous sections pro- 
duces 19 sets of bandpowers and covariance matrices, one 
from each field. In either the limit of equal noise or the 
limit of sample variance domination, the optimal weight 
for each field would be its effective area (i.e., the inte- 
gral of its window). Since we are approximately in these 
limits, we use area-based weights. Thus the combined 
bandpowers and covariance matrix are given by 



Dt, 



w 



where 



.4' 



w 



(10) 



(ii) 



(12) 



is the area-based weight of the i field. The area A 1 is 
the sum of the window for the i th field. 



We calculate the final covariance matrix as the sum of 
the signal plus noise covariance matrix, the beam covari- 
ance matrix, and the calibration covariance matrix. For 
the signal and noise terms, we combine the signal plus 
noise covariance matrices from all fields using Equation 
|TTj We condit ion this combined covariance matrix using 
Equation 11 from Kll. For the beam covariance matrix, 
we t ake the beam covariance matrices for each year (see 
^3.5), and combine them into one composite beam co- 



variance matrix using the area-based weight scheme. In 
this step, we take care to account for the beam errors 
that are correlated between years. Finally, we add the 
calibration covariance matrix, defined as C™' = e 2 DiDj, 
where e = 0.027 is the 2.6 % u ncertainty in the SPT power 
calibration discussed in 92.4 



3.7. Bandpower Window Functions 

Bandpower window functions are necessary to compare 
the measured bandpowers to a theoretical power spec- 
trum. The window function Wf II is defined as 



th 



(13) 



Following the formalism described in Section |3.4[ we 
can write this as 



PveMtuFtBiC, 



)W r Vt 



which implies that 

W\/l = {K-^wPvt'MetFtBe. 



th 



(14) 



(15) 



We calculate the bandpower window functions to be 
used for the final spectrum measurement as the weighted 
average of the bandpower window functions from each 
field. 

4. TESTS FOR SYSTEMATIC ERRORS 

It is important to verify that the data are unbiased 
by systematic errors. As is common in CMB analyses, 
we check for possible systematic errors by performing a 
suite of null tests which are frequently referred to as jack- 
knife tests. In each null test, all observations are divided 
into two equally sized sets based on a possible source of 
systematic error. Difference maps are then calculated by 
subtracting the two sets, thus removing the astrophysical 
signal. The power spectrum of the difference maps is cal- 
culated as described in the last section. This spectrum 
is compared to an "expectation spectrum," the power 
we expect to see in the absence of contamination from 
systematic errors. The expectation spectrum will gener- 
ically be non-zero due to small differences in observation 
weights, filtering, etc, and is calculated by applying the 
null test to simulated maps. The expected power is small 
{Di < 2/iK 2 at all multipoles) for all tests. 

We perform six null tests: 

• Time: Observations are ordered by time, then di- 
vided into first-half and second-half sets. This 
tests for long-term temporally varying systematic 
effects. 

• Scan Direction: Observations are divided into 
maps made from left-going scans and right-going 
scans. This tests for scan-synchronous and scan- 
direction-dependent systematic errors. 
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• Azimuthal Range: We split the data into ob- 
servations taken at azimuths that we expect to 
be more or less susceptible to ground pickup. 
These azimuth ranges are determined from maps 
of the 2009 data made using "ground-centered" 
(Azimuth/Elevation) coordinates in which ground 
pickup adds coherently, as opposed to the usual 
"sky-centered" (R.A./decl.) coordinates. We use 
the ground-centered maps that were made for the 
analysis presented by Kll. Although we detect 
emission from the ground on large scales {£ ~ 50) 
in these ground-centered maps, this is not expected 
to bias our measurement; the amplitude of the 
ground pickup is significantly lower on the smaller 
angular scales for which the bandpowers are be- 
ing reported, and the observations for a given field 
are distributed randomly in azimuth. We use the 
azimuth-based null test to verify this assertion. 

• Moon: Observations are divided into groups based 
on when the Moon was above and below the hori- 
zon. 

• Sun: Observations are divided into groups based 
on when the Sun was above and below the horizon. 
In this test, we only include fields in which more 
than 25% of the observations were taken with the 
Sun above the horizon. 

• Summed Bolometer Weights: We calculate the sum 
of all bolometer weights during each observation 
and order maps based on this sum. This tests for 
bias introduced by incorrectly weighting observa- 
tions or incomplete coverage in some maps. 

For each test, the x 2 of the residual power is calcu- 
lated relative to the expectation spectrum in five bins 
with S£ = 500. We calculate the probability to exceed 
(PTE) this value of x 2 for five degrees of freedom. All 
null tests had reasonable PTE's, as listed below, with 
the exception of the Azimuthal Range null test, which 
produced a low PTE for the original set of observations. 
This was interpreted as evidence for some ground con- 
tamination. This interpretation was tested by cutting 
several sets of 5% of the data, and re-calculating the 
Azimuthal Range null test. Removing random 5% sets 
of the data did not change the failure of the Azimuthal 
Range null test. However, cutting the 5% of the obser- 
vations from each field with the highest expected ground 
contamination resulted in passing the Azimuthal Range 
null test. Thus this cut was incl uded with the other 
observation cuts, as described in S 12.11 It is worth not- 



ing that the Azimuthal Range null test is the worst-case 
scenario for ground pickup; this null test systematically 
aligns azimuth ranges to maximize the ground contam- 
ination. In the analysis of the power spectrum of the 
sky signal, the ground signal will add incoherently as the 
azimuth changes, thus reducing the power from ground 
contamination to a much lower level than in this null 
test. 

We find a flat distribution of PTEs ranging from 0.02 
to 0.99 for individual fields. The combined PTEs for the 
Time, Scan Direction, Azimuthal Range, Moon, Sun, and 
Summed Bolometer Weights are 0.26, 0.14, 0.17, 0.13, 
0.30, 0.63, respectively. It is important to note that 



these null tests are extremely conservative for the SPT 
power spectrum where the uncertainties are sample vari- 
ance dominated over most of the range reported in this 
work. It is possible to have a failure in these null tests 
without a significant impact on the final power spectrum. 
The measured power Dg in each null test was less than 
2.2 fiK 2 in all bins. We thus find no significant evidence 
for systematic contamination of SPT bandpowers. 

5. BANDPOWERS 

Following the analysis presented in [|3j we measure the 
CMB temperature anisotropy power spectrum from 2540 
deg 2 of sky observed by the SPT between 2008 and 2011. 
We report bandpowers in bins of S£ — 50 between 650 < 
I < 3000. The bandpowers and associated errors are 
listed in Table [2] and shown in Figures [3]and[4| The 
bandpowers, covariance matrix, and window functions 
are available for download on the SPT website. 28 

These bandpowers clearly show the third to ninth 
acoustic peaks. As Figure|4]demonstrates, the anisotropy 
power measured by this analysis (at 150 GHz, with a 
50mJy point-source cut) is dominated by primary CMB, 
with secondary anisotropy and foregrounds contributing 
significantly only at the highest multipoles. These band- 
powers provide the most precise measurement to date of 
the CMB power spectrum over the entire multipole range 
presented in this analysis. 

6. COSMOLOGICAL CONSTRAINTS 

The SPT bandpowers are high signal-to-noise measure- 
ments of the CMB temperature anisotropy over a large 
range of angular scales, and can be used to perform sen- 
sitive tests of cosmological models. In this section, we 
present the constraints these bandpowers place on cosmo- 
logical models. We first constrain the standard ACDM 
cosmological model. Next, we extend this model to con- 
strain the amplitude of gravitational lensing of the CMB. 
We then consider models with free spatial curvature and 
constrain the mean curvature of the observable universe. 
We finally consider tensor perturbations and discuss the 
implications of our observations for simple models of in- 
flation. A wider range of cosmological models are tested 
in a companion paper H12. 

We parametrize the ACDM model with six parame- 
ters: the baryon density fib/i 2 , the density of cold dark 
matter fl c h 2 , the optical depth to reionization r, the an- 
gular scale of the sound horizon at last scattering 9 S , the 
amplitude of the primordial scalar fluctuations (at pivot 
scale ko = 0.05 Mpc -1 ) A^, and the spectral inde x of 
the scalar fluctuations n s . With the exception of § |6.8[ 
we consider only flat-universe models where the mean 
curvature of the observable universe VL^ = 0. In addition 
to the six parameters described above, we report sev- 
eral derived parameters that are calculated from the six 
ACDM parameters. These are the dark energy density 
Q\, the Hubble constant Hq in units of km s _1 Mpc -1 , 
the current amplitude of linear matter fluctuations eg 
on scales of 8ft" 1 Mpc, the redshift of matter- radiation 
equality zeq, and a hybrid distance ratio reported by 
baryon acoustic oscillation (BAO) experiments at two 
different redshifts r s /D v (z = 0.35) and r s /D v (z — 0.57), 
where r s is the comoving sound horizon size at the baryon 

28 http:/ /pole. uchicago.edu/public/data/storyl2/ 
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TABLE 2 

SPT Bandpowers and Bandpower Errors 



I range 


4ff 


Di [fiK 2 ] 


a [ M K 2 ] 


£ range 


4a 


D e [ M K 2 ] 


a [/xK 2 ] 


651 


- 700 


671 


1786.2 


59.5 


1851 - 


1900 


1865 


276.1 


5.2 


701 


- 750 


720 


1939.3 


66.9 


1901 - 


1950 


1915 


238.2 


4.5 


751 


- 800 


770 


2426.4 


67.2 


1951 - 


2000 


1966 


242.8 


4.6 


801 


- 850 


820 


2577.1 


68.3 


2001 - 


2050 


2015 


245.8 


4.9 


851 


- 900 


870 


2162.3 


53.8 


2051 - 


2100 


2064 


229.9 


4.5 


901 


- 950 


920 


1588.8 


39.0 


2101 - 


2150 


2114 


194.2 


3.8 


951 


- 1000 


969 


1144.3 


29.6 


2151 - 


2200 


2164 


170.6 


3.4 


1001 


- 1050 


1019 


1068.0 


27.2 


2201 - 


2250 


2213 


140.5 


2.8 


1051 


- 1100 


1069 


1215.7 


28.5 


2251 - 


2300 


2265 


135.0 


2.6 


1101 


- 1150 


1118 


1193.8 


29.1 


2301 - 


2350 


2313 


128.3 


2.4 


1151 


- 1200 


1169 


1141.1 


29.8 


2351 - 


2400 


2364 


124.8 


2.7 


1201 


- 1250 


1218 


924.8 


23.1 


2401 - 


2450 


2413 


115.8 


2.2 


1251 


- 1300 


1269 


771.7 


17.9 


2451 - 


2500 


2462 


100.7 


2.2 


1301 


- 1350 


1318 


723.1 


17.7 


2501 - 


2550 


2512 


96.7 


2.3 


1351 


- 1400 


1367 


754.6 


16.6 


2551 - 


2600 


2562 


83.3 


2.0 


1401 


- 1450 


1417 


847.3 


17.0 


2601 - 


2650 


2613 


85.7 


1.8 


1451 


- 1500 


1468 


718.7 


13.8 


2651 - 


2700 


2663 


83.9 


1.9 


1501 


- 1550 


1517 


625.0 


11.3 


2701 - 


2750 


2712 


76.4 


1.8 


1551 


- 1600 


1567 


468.1 


10.2 


2751 - 


2800 


2761 


71.7 


1.8 


1601 


- 1650 


1617 


395.7 


7.9 


2801 - 


2850 


2811 


62.9 


1.7 


1651 


- 1700 


1666 


390.7 


7.0 


2851 - 


2900 


2860 


57.6 


1.6 


1701 


- 1750 


1717 


396.6 


6.9 


2901 - 


2950 


2910 


57.6 


1.6 


1751 


- 1800 


1766 


390.7 


6.9 


2951 - 


3000 


2961 


56.6 


1.6 


1801 


- 1850 


1815 


336.7 


6.2 













Note. — The ^-band range, weighted multipole value l e g, bandpower De, and associated bandpower uncertainty a of the 
SPT power spectrum. The errors are the square-root of the diagonal elements of the covariance matrix, and do not include 
beam or calibration uncertainties. 



drag epoch, D v (z) = [(1 + z) 2 D%{z)cz/H(z)} 1 ' 3 , D A (z) 
is the angular diameter distance, and H(z) is the Hubble 
parameter. 

6.1. Foreground treatment 

We marginalize over three foreground terms in all pa- 
rameter fitting. We have tested that all parameter con- 
straints are insensitive to the details of the assumed fore- 
ground priors. The total foreground power, , can be 
expressed as: 



D. 



D'i 



+ Df 

These two terms represent the following: 



(16) 



• Power from galaxies (Df a ), which can be subdi- 
vided into a clustering term and a Poisson term. 
The Gaussian priors used by Kll are applied to 
the amplitude of each term at I = 3000. For the 
clustering term, the prior is D^nnn = 5.0 ± 25 fjK 



(2011) 



based on measurements by Shirokoff et al. 
The angular dependence or the clustering term is 
Df h oc £°' s , which has been modified from that 
assumed by Kll to agree better with r ecent mea- 
surements (e.g., |Addison et al. (2012 1; |Reichardt 
For the Poisson term 

2 



et al. 



(2012b)). For the Poisson term, the prior 
, (ill0 — 19.3 ± 3.5^iK 2 . This is based on the 
power from sources with 5 i 5oqhz < 6.4 mJy, as 



measured in Shirokoff et al. (20111, and the power 



as measured in Vieira et al. (2010) and Marriage 
et al. (2011). The Poisson term is constant in C'e 



and thus varies as D^ s oc £ 2 . 

SZ power (Df z ). The thermal and kinetic SZ ef- 
fects are expected to contribute to the observed 
CMB temperature anisotropy. Both effects are ex- 
pected to have similar power spectrum shapes over 
the angular scales relevant to this analysis. There- 
fore we adopt a single template to describe both ef- 
fects. The chosen template is the thermal SZ model 
We set a Gaussian prior on 



from Shawetal. (2010 



the amplitude of this term of Z)fp m = 5.5±3.0 [iK z 



as measured in Shirokoff et al. (2011). This ampli 



from sources with 6.4 mJy < SxsoGHz < 50 mJy, 



tude is defined at 153 GHz, corresponding to the 
effective SPT band center. 

In the above model, we have not accounted for the 
emission from cirrus-like dust clouds in the Milky Way. 
Repeating the calculation performed in Kll, we cross- 
correlate the SPT maps with predictions for the galactic 
d ust emission at 150 GHz in the SPT fields using model 8 
of Finkbeiner et al. (1999). We use this cross-correlation 
to estimate the power from the galactic dust in the SPT 
fields, and find that it is small compared to the primary 
CMB power and the SPT bandpower errors. Specifically, 
subtracting this cirrus power can be balanced by moving 
the foreground terms by amounts that are small com- 
pared to their priors, so that the change in \ 2 is l ess 
than 0.1 o. Thus we conclude that galactic dust does not 
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Fig. 3. — Left panel: The SPT power spectrum. The leftmost peak at I ~ 800 is the third acoustic peak. Right panel: A compa rison 
of the new S PT bandpowers with other recent measurements of the CMB damping tail from ACBAR (Reichardt ct al. 2009), ACT | |Das| 
|et al,|2011b| , and SPT (Kll). Note that the point source masking threshold differs between these experiments which can attect the power 
at the highest multipoles. In order to highlight the acoustic peak structure of the damping tail, we plot the bandpowers in the right panel 
as t 4 Ce/(2w), as opposed to Di = £(£ + l)C^/(27r) in the left panel. The solid line shows the theory spectrum for the ACDM model + 
foregrounds that provides the best fit to the SPT+ WMAP7 data. The bandpower errors shown in these plots contain sample and noise 
variance terms only; they do not include beam or calibration uncertainties. 



significantly contaminate the SPT power spectrum. 

6.2. Estimating Cosmological Parameters 

Our baseline model contains nine parameters: six for 
the primary CMB and three for foregrounds. We explore 
the nine-dimensional parameter space usin g a Markov 
Chain Mon te Carlo (MCMC) technique ( Christensen 
et al.||2001[) implemented in the CosmoMC^ (Lewis & 
UridlepUU^I ) software package. For reason s of speed, we 
use PICO 3 " dFendt fc Wandeltl |2007b|a|), trained with 
CAMB 31 ( [Lewis et al.|2000[ ), to calculate the CMB power 
spectrum. We have trained PICO for a ten-parameter 
model that includes ACDM as well as several exten- 
sions. We use PICO when working with any subset of 
this model space, and CAMB for all other extensions to 
ACDM. The effects of gravitational lensing on the power 
spectrum of the CMB are cal culated using a cosmology - 
dependent lensing potential ( Lewis fc Challinor||2006 [). 
To sample the posterior probability distribution m re- 
gions of very low probability, we run "high-temperature" 
chains in which the true posterior, P, is replaced in the 
Metropolis Hastings algorithm by Pt — P 1 ^ 6 . This 
allows the chain to sample the parameter space more 
broadly. We recover the correct posterior from the chain 
by importance sampling each sample with weight Pj Pt- 

6.3. Goodness of fit to the ACDM Model 

We quantify the goodness of fit of the ACDM model to 
the SPT bandpowers by finding the spectrum which best 
fits the SPT bandpowers and calculating the reduced x 2 
for the SPT data. The reduced x 2 for the SPT data is 
45.9/39 (PTE=0.21), thus the ACDM model is a good 
fit to the SPT bandpowers. In H12, we consider several 



extensions to the ACDM model, and find that the data 
show some preference for several of those extensions. 

6.4. External Datasets 

In this work, we focus on parameter constraints from 
the CMB data, sometimes in conjunction with measure- 
ments of the Hubble constant (Ho) or the BAO feature. 
For CMB measurements, we use the SPT bandpowers 
presented here as well as the WMAP bandpowers pre- 
sented in WMAP7. For Hq m easurements, we us e the 
low-redshift measurement from jRiess et al. (20111. For 
the BAO feature, we use a combination of three measure- 
ments at different redshifts: the WiggleZ survey cov ering 



the redshift range 0.3 < z < 0.9 (|Blake et al.pbllll .Jhe 
SDSS-II survey (DR7) covering 0.16 < z < 0.44 l|Pad 



manabhan et a l.l|2012D, and the BO SS survey covering 



29 http : / /cosmologist . info / cosmomc / 

30 https : //sites . google . com/a/ucdavis . edu/pico 

31 



http://camb.info/ (January 2012 version) 



0.43 < z < 0.7 ( |Anderson et al.||2012[ ). 

Before combining the CMB, Ho, and BAO datasets, we 
check their relative consistency within the ACDM model. 
We quantify this consistency by calculating the x 2 n i n us ~ 
ing a reference dataset (e.g., CMB) and comparing it to 
the Xmin obtained using a new dataset (e.g., CMB+H ). 
For example, xL n ,[CMB+H ] ~ Xmin,[CMB]= - 08 - The 
probability to exceed this Ax 2 given the one new degree 
of freedom provided by the Hq measurement is 0.78, cor- 
responding to an effective Gaussian significance of 0.3 ct. 
Using this metric, we find that 

• CMB and H differ by 0.3 ct. 

• CMB and BAO differ by 1.5 ct. 

• (CMB+BAO) and H differ by 1.8 ct. 

• (CMB+H ) and BAO differ by 2.1 ct. 

There is some tension between these datasets in the 
context of the ACDM model. This could be evidence 
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Fig. 4. — The SPT bandpowers (blue), WMAP7 bandpowers (orange), and the lensed ACDM+foregrounds theory spectrum that 
provides the best fit to the SPT+ WMAP7 data shown for the CMB-only component (dashed line), and the CMB+foregrounds spectrum 
(solid line). As in Figure [3] the bandpower errors shown in this plot do not include beam or calibration uncertainties. 

all parameters; 9 S changes the most significantly among 
the five free ACDM parameters, moving by 1.5 a and 
tightening by a factor of 2.2 from WMAP7 to SPT. The 
SPT bandpowers measure 6 S extremely well by virtue of 
the sheer number of acoustic peaks - seven - measured 
by the SPT bandpowers. The SPT constraint on n s is 
broader than the constraint from WMAP7 due to the 
fact that WMAP7 probes a much greater dynamic range 
of angular scales. Degeneracies with n s degrade the SPT 
constraints on A%, the baryon density and, to a lesser 
extent, the dark matter density. 

6.6. Combined ACDM constraints 

Next, we present the constraints on the ACDM 
model from the combination of SPT and WMAP7 data. 
As previously mentioned, we will refer to the joint 
SPT+ WMAP7 likelihood as the CMB likelihood.' We 
then extend the discussion to include constraints from 
CMB data in combination with BAO and/or Hq data. 

We present the CMB constraints on the six ACDM 
parameters in the fourth column of Table [3] Adding 



for a departure from ACDM, a systematic error in one 
or more of the data sets, or simply a statistical fluctua- 
tion. We assume the uncertainties reported for each of 
the datasets are correct and combine them to produce 
many of the results presented here. 

6.5. SPT-only ACDM constraints 

We begin by examining parameter constraints from the 
SPT bandpowers alone. The SPT-only parameter con- 
straints provide an independent test of ACDM cosmology 
and allow for consistency checks between the SPT data 
and other datasets. Because the scalar amplitude A|, 
and the optical depth r are completely degenerate for 
the SPT bandpowers, we impose a WMAP7-based prior 
of r = 0.088 ± 0.015 for the SPT-only constraints. 

We present the constraints on the ACDM model from 
SPT and WMAP7 data in columns two to four of Table0 
As shown in Figure pi the SPT bandpowers (including 
a prior on r from WMAP7) constrain the ACDM pa- 
rameters approximately as well as WMAP7. The SPT 
and WMAP7 parameter constraints are consistent for 
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Fig. 5. — The one-dimensional marginalized likelihoods of the six parameters of the ACDM model, plus two derived parameters: the dark 
energy density and the Hubble constant Ho- The constraints are shown for the SPT-only (blue dot-dashed lines), WMAP7-oialy 
(red dashed lines), and SPT+ WMAP7 (black solid lines) datasets. With the exception of t, the SPT bandpowers constrain the 
parameters approximately as well as WMAP7 alone. In particular, the SPT bandpowers measure the angular sound horizon 8 a extremely 
well because they measure seven acoustic peaks. In the SPT-only constraints, the WMAP7 measurement of r has been applied as a prior; 
because of this we do not plot an SPT-only line on the r plot. 



TABLE 3 

ACDM Parameter Constraints from the CMB and external datasets 



Parameter 


WMAP7 


SPT< a ) 


CMB 
(SPT+ WMAP7) 


CMB+iio 


CMB+BAO 


CMB+H ±BAO 


Baseline paramete 
100 U b h 2 
U c h 2 
10 9 A| 

n s 

ioo e s 

T 


•s 

2.231 ± 0.055 
0.1128 ±0.0056 

2.197 ±0.077 

0.967 ±0.014 
1.0396 ±0.0027 

0.087 ±0.015 


2.30 ±0.11 
0.1056 ±0.0072 
2.164 ± 0.097 
0.926 ±0.029 
1.0441 ± 0.0012 
0.087 ± 0.015 


2.229 ±0.037 
0.1093 ±0.0040 

2.142 ±0.061 
0.9623 ± 0.0097 
1.0429 ±0.0010 

0.083 ±0.014 


2.233 ±0.035 
0.1083 ±0.0033 

2.138 ±0.062 
0.9638 ± 0.0090 
1.0430 ± 0.0010 

0.084 ± 0.014 


2.204 ± 0.034 
0.1169 ±0.0020 

2.161 ± 0.057 
0.9515 ±0.0082 
1.04215 ± 0.00098 

0.076 ±0.012 


2.214 ±0.034 
0.1159 ±0.0019 

2.160 ±0.057 
0.9538 ± 0.0081 
1.04236 ± 0.00097 

0.077 ±0.013 


Derived parameter 

n A 

H 

^EQ 

100-jjMz = 0.35) 
100^(z = 0.57) 


0.724 ±0.029 
70.0 ± 2.4 

0.819 ±0.031 
3230 ± 130 
11.43 ±0.37 
7.58 ±0.21 


0.772 ±0.033 
75.0 ±3.5 

0.772 ±0.035 
3080 ± 170 
12.15 ±0.55 
7.98 ±0.31 


0.750 ±0.020 
72.5 ± 1.9 

0.795 ±0.022 
3146 ± 95 

11.81 ± 0.29 
7.80 ±0.16 


0.755 ±0.016 
73.0 ± 1.5 

0.791 ± 0.019 
3124 ± 78 

11.89 ±0.24 
7.84 ± 0.13 


0.709 ±0.011 
69.11 ± 0.85 

0.827 ± 0.015 

3323 ± 50 
11.28 ±0.12 

7.505 ±0.068 


0.7152 ± 0.0098 
69.62 ± 0.79 
0.823 ±0.015 
3301 ± 47 
11.35 ± 0.12 
7.545 ± 0.065 



Notes: The constraints on cosmological parameters from the ACDM model, given five different combinations of datasets. 
We report the median of the likelihood distribution and the symmetric 68.3% confidence interval about the mean. 



(a) We impose a IT r M,4P7-based prior of r = 0.088 ± 0.015 for the SPT-only constraints. 

^ Derived parameters are calculated from the baseline parameters in CosmoMC. They are defined at the end of § 6.1 
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SPT bandpowers to the WMAP7 data tightens these 
parameter constraints considerably relative to WMAP7 
alone. Of these parameters, the constraint on 8 S sees the 
largest improvement; adding SPT data decreases the un- 
certainty on 6 S by a factor of 2.70 relative to WMAP7 
alone. Constraints on on 0;,/i 2 , f2 c /i 2 , and fl\ tighten by 
factors of 1.49, 1.40, and 1.45, respectively. For compari- 
son, the addition of the Kll bandpowers to WMAP7 led 
to improvements of 1.33, 1.16, and 1.11, respectively. Fi- 
nally, the constraint on the scalar spectral index tightens 
by a factor of 1.44 to give n s < 1.0 at 3.9 a. 

The preferred values for fl b h 2 and fL c h 2 for WMAP7 
do not shift significantly with the addition of the SPT 
data and, therefore, neither does the sound horizon, r s , 
which depends only on these parameters in the ACDM 
model. Thus the shift in 9 S = r s /DA, driven by the 
SPT acoustic peak locations, must lead to a shift in Da- 
Shifting Da requires shifting Q\ (or, equivalently, H ). 
This shift in J7a and Ho can be seen in Figure [5 

We show the parameter constraints after adding the 
Hq and/or BAO data to the CMB data in the last three 
columns of Table [3] Combining the CMB bandpow- 
ers with this low-redshift information tightens the con- 
straints on Q, c h 2 and fl\ by a further factor of 1.2 or 2 for 
CMB+Ho and CMB+BAO, respectively, with smaller 
but significant improvements to other parameters. Of 
special note is the constraint on the scalar spectral index 
which tightens to n s = 0.9638 ±0.0090 for the CMB+H 
dataset, 0.9515±0.0082 for the CMB+BAO dataset, and 
0.9538 ± 0.0081 for the CMB+£T +BAO dataset. These 
constraints correspond to a preference for n s < 1 at 4.0 a, 
6.1 cr, and 5.9 a respectively, for these three data combi- 
nations. This is the most signi fican t reported measure- 
ment of n s 7^ 1 to date. See § |6.9| for a more detailed 
discussion of constraints on n s . 

6.7. Gravitational Lensing 

As CMB photons travel from the surface of last scat- 
tering to the Earth, their paths are deflected by gravita- 
tional interactions with intervening matter. This grav- 
itational lensing encodes information about the distri- 
bution of matter along the line of sight, providing a 
probe of the distance scale and growth of structure at 
intermediate redshifts (0.5 < z < 4). Lensing distorts 
the CMB anisotropy by shifting the apparent position 
of CMB photons on the sky, with typical deflection an- 
gles of 2.5 arcminutes which are coherent over degree 
scales. This process mixes power between multipoles in 
the CMB temperature power spectrum, which smooths 
the acoustic peak structure and increases the power in 
the damping ta il at small angular scales (see Lewis & 
Chall inor||2006 for a review) . 

The ACDM model already includes the effects of grav- 
itational lensing. To quantify the sensitivity of the 
SPT bandpowers to gravitational lensing, we extend the 
ACDM model to include one additional free parameter, 
Al, which rescales the lensing potential power spectrum, 
, according to 



a 



ArCt 



(17) 



We re-calculate Cf in a cosmology-dependent manner 
at each point in the MCMC. In effect, the Al parameter 
modulates the amplitude of gravitational lensing. Set- 
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Fig. 6. — The SPT bandpowers allow a significant detection of 
gravitational lensing through the effective smoothing of the acous- 
tic peaks. Here, we show the one-dimensional likelihood function 
for Al , a rescaling parameter for the gravitational lensing potential 
power spectrum (Cf -> A L Cf' p ). The SPT+ WMAP7 data lead 
to a 8.1 ff detection of CMB lensing, the most significant detection 
to date. 

ting Al — 1 corresponds to the standard theoretical pre- 
diction and recovers the standard ACDM model, while 
setting Al = corresponds to no gravitational lensing. 
In the parameter fits, the range of Al is allowed to ex- 
tend well above 1 and below 0. 
The first detection of gravitational len sing in the 



CMB used lensing-galaxy cross-correlations 
20071, and subsequent papers using this tec 



achieved higher signal-to- noise detections (Bleem et al 



2012 Sherwin et al. 2012 ) 



Smith et al 
mique have 



"The impact ot lensing on the CMB power spectrum has 
been detected in combinations of WMAP with ACBAR 



iReichardt et al. 2009J), WMAP with ACT ( Das et al . 
blla|, and WM APwth SPT (Kll). Using this effect 



Das et al. 



1.3±g J at 68% confidence. 



2011a[ ) found A L 
Kll found that the constraint on A - 65 had the most 
mssian shape and thus reported A L 65 = 0.94 ± 0.15, a 
5 a detection of lensing. 

CMB lensing has also been detected throu gh the CMB 



tempe rature four- point function in ACT ( Das et al 
201 lbH and SP T ( jvan Engelen et al.|[2012l ) data. JDas 

^"TTT dom 1 Kh „ m J APT +^ m „ nm ,v n 32 ^ME ~ 



(2011M used ACT data to measure 
0.29. Inlvan Engelen et all {l2012h, the S 



et al.| 

1.16 ± 0.29. lh |van Engelen et al.N2012[ ), the SPT four- 
point analysis was applied to a subset of the data used 
in this work to measure Al — 0.90 ± 0.19, which was 
previously the most significant detection of CMB lensing 
to date, ruling out no lensing at 6.3 a. 
We determine the significance of the observed CMB 



32 In |Das et al.|{2011b| , Al is calculated as the best-fit am- 
plitude to the lensing potential in the cosmological model with 
the maximum likelihood ACDM parameters. This is in contrast to 
what was done for the CMB temperature power spectrum measure- 
ments of Al , where constraints on Al have been marginalized over 
cosmological parame ters. The correspo nding maximum likelihood 
measure from | van Engelen et a~ | 2012| is A\^ = 0.86 ± 0.16. 



14 



lensing by constraining Al with the measured CMB 
power spectrum. The significance of the detection is 
quantified by calculating the probability for Al < 0, 
P(A L < 0). As A h — is far out in the tail of the 
likelihood distribution of this parameter, we use high- 
temperature MCMC's chains to estimate P(Al < 0). 
Using SPT data only, we measure P^ SPT ^(A L < 0) < 
1.3 x 10 -9 , the equivalent of a 5.9 ct preference for Al > 
in a Gaussian distribution. For SPT+ WMAP7 we mea- 
sure 

P {CMB} (A L < 0) < 2.4 x 1(T 16 , (18) 

which corresponds to a 8.1 ct detection of lensing in a 
Gaussian distribution. 

Next, we report constraints on Al- Using 
SPT+ WMAP7 data we find 



A, 



0.86 



+0.15(+0.30) 
-0.13(-0.25) 



(19) 



where asymmetric 1 ct (68.3%) and 2 ct (95.5%) errors are 
shown. The observed lensing amplitude is consistent at 
1 ct with theoretical predictions in the ACDM model. 

6.8. Mean Curvature of the Observable universe 

The low-redshift information imprinted on the CMB 
by gravitational lensing, along with the other informa- 
tion in the CMB anisotropy power spectrum, enables the 
placement of tight constraints on the mean curvature of 
the observable Universe. The magnitude of the mean 
curvature today can be parameterized by Qk = —K/Hq 
where -y/l/l-K"! is the length scale over which departures 
from Euclidean geometry become important. Infl ation- 
ary m odels generically predict |$lfc| < 10 -5 (e.g., Knox 



2006 1 ; thus a significant measurement of flk ^ would 
challenge our standard picture of the very early universe. 
A positive value for Slk(K < 0) could be obtained by the 
nucleation of a bubble of lower vacuum energy in a sur- 
rounding medium with higher vacuum energy follo wed 
by a short period of inflation (Bucher et al. 19951. A 
determination that Sl k is negative with high statistical 
significance would be very interesting; such a detection 
would be difficult to understand in the theoretical frame- 
work of inflation, challenge the string theory landscape 
pi cture, and rule o ut the de Sitter equilibrium cosmology 
of |Albrecht| ( |2011| . 

Absent lensing effects, one can leave the CMB power 
spectrum unchanged while simultaneously varying flk 
and f?A in a way that keeps the distance to last scat- 
tering fixed ( |Bond et al.|l997||Zaldarriaga et al.|l997j ). 33 
Historically the CMB data placed very coarse constraints 
on fifc, with finer constraints only possible with the ad- 
dition of other data sensitive to Slf. and JIa s uch as Hp 



measurements and determinations of Sl m (e.g., Dodelson 
fe Knox||2000 ) from, for example, the baryon fraction in 
clusters of galaxies ( White et al.||1993 ). 

The sensitivity of the CMB to low-redshift informa- 
tion through gravitational lensing makes it possible to 
constrain the mean curvature of the observable universe, 
and thus the cosmological constant, using the CMB 
alone. The lensing amplitude is sensitive to the dis- 
tance and growth of structure at intermediate redshifts 

33 There is an exception to this at very large scales due to the late 
ISW effect, but sample variance makes these changes unobservably 
small in the CMB power spectrum. 



(0.5 < z < 4). These observables are, in turn, sensi- 
tive to curvature, dark energy, and neutrino masses, as 
discussed in H12. The recent detections of CMB lens- 
ing have measured a n amplitude that is consistent with 
£7fi-~0 and S1a~0.7 flDas et al.|p011a|b[ Shcrwin et al. 
2011||Keisler et al.||20Tl[ jvan Engelen et aT^p[Jl2] . Sim- 
ply put, the strength of CMB lensing in a universe with 
no dark energy and positive mean curva ture would be 
much larger than that observed (see e.g., Sherwin et al. 
20lT| . 

Using the SPT+ WMAP7 bandpowers, we measure the 
mean curvature of the observable universe using only the 
CMB: 



-0.003 



+0.014 
-0.018 ' 



(20) 



This result tightens curvature constraints over WMAP7 
combined with low-redshift probes by ~20%. This con- 
straint is consistent with zero mean curvature, and cor- 
responds to a dark energy density of JIa = 0.740l ; 054 
and Hubble constant, H Q = 70.9±g o (km s" 1 Mpc" 1 ). 
This measurement rules out Oa = at 5.4 ct using the 
CMB alone. The right panel of Figure [7] shows the cor- 
responding two-dimensional marginalized constraints on 
SIm and SI a- We have confirmed that the strength of 
this constraint relies on the lensing signal; allowing Al 
to vary freely causes the curvature constraint to degrade 
dramatically. 

CMB lensing enables an independent constraint on 
curvature, although the most powerful curvature con- 
straints still come from combining CMB data with other 
low-redshift probes (e.g., Hq, BAO). The curvature con- 
straint using CMB+# data is Q k = 0.0018 ± 0.0048, 
while the constraint using CMB+BAO data is Sl k = 
-0.0089 ± 0.0043. The tightest constraint on the mean 
curvature that we consider comes from combining the 
CMB, H , and BAO datasets: 



fi« 



-0.0059 ± 0.0040 . 



(21) 



While the CMB+BAO constraint shows a 2.0 ct pref- 
erence for ttj. < 0, the significance of this preference de- 
creases as more data are added. The tightest constraint, 
coming from CMB+ff +BAO, is consistent with zero 
mean curvature at 1.5 ct. These results are summarized 
in Figure [7J As discussed by H12, other extensions of 
ACDM can also explain the data (e.g., allowing for both 
non-zero mean curvature and a dark energy equation of 
state w =/= —1), thus these constraints are significantly 
degraded when multiple extensions to ACDM are simul- 
taneously considered. 

6.9. Inflation 

Cosmic infla tion is an accelerated expansion in the 
early univer se ( |Guth|1981| |Linde|1982| |Albrecht fc Stein 



hardt 19821 that generically leads to a universe with 



nearly zero mean curvature and a nearly sca le-invariant 
spectrum of "initial" density perturbations jM ukhanov 
fc Chibisov(l98Tl |Hawking|1982| |Starobinsky|lO^| |Gut,n 



Pi||1982J [Bardeen et al.||1983p that evolved to pro 



duce the observed spectrum of CMB anisotropies. Mod- 
els of inflation compatible with current data generally 
predict, over the range of observable scales, scalar and 
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Fig. 7. — The low-redshift information imprinted on the CMB by gravitational lensing, along with the other information in the CMB 
anisotropy power spectrum, enables the placement of tight constraints on the mean curvature of the observable Universe. The addition of 
low-redshift probes further tighten CMB-only constraints on the mean curvature. Left panel: The one-dimensional marginalized constraints 
on Cl k from SPT+ WMAP7 (black solid line), SPT+ WMAP7+Ho (orange dashed line), SPT+ WMAP7+BAO (black dotted line), 
and SPT+ WMAP7+BAO+H (blue dot-dashed line). The SPT+ WMAP7 datasets measure the mean curvature of the observable 
universe to a precision of ~1.5%, while combining SPT, WMAP7, and either Ho or BAO data reduces the uncertainty by a factor of ~3. 
Right panel: The two-dimensional constraints on Qm and from the SPT+ WMAP7 data alone. The SPT+ WMAP7 data rule out 
Q A = at 5.4 a. 



TABLE 4 

Constraints on n s and r from CMB and external datasets 



Model 


Parameter 


CMB 


CMB±ii"o 


CMB±BAO 


CMB±# ±BAO 






(SPT+ WMAP7) 








ACDM 


n s 


0.9623 ±0.0097 


0.9638 ± 0.0090 


0.9515 ± 0.0082 


0.9538 ± 0.0081 


ACDM+r 


n s 


0.969 ±0.011 


0.9702 ± 0.0097 


0.9553 ± 0.0084 


0.9577 ±0.0084 




r (95% C.L.) 


< 0.18 


< 0.18 


< 0.11 


< 0.11 



tensor perturbations well characterized by a power law 
in wavenumber fc, 



Ai(fc) = A 2 R (k Q ) 



Al(k) = A 2 h (k ) 



ko 



n.-i 



n, 



(22) 



(23) 



Here A^(fco) is the amplitude of scalar (density) pertur- 
bations specified at pivot scale k — fco, with scale de- 
pendence controlled by the index n s , while A|(fco) is the 
amplitude of tensor (gravitational wave) perturbations 
specified at the same pivot scale, with scale dependence 
set by nt ■ The amplitude of the tensor perturbation spec- 
trum is expressed in terms of the tensor-to-scalar ratio 



A^(fc) 



A|(fc) 



(24) 



fc=0.002 Mpc- 



For single-field models in slow-roll i nflation, n« and r are 
related by a consistency equation ( |Copeland et al.|[l993 



Kinney et aL||2 008): 



-r/8. 



(25) 



The tensor and scalar perturbations predicted by such 
models of inflation can thus be characterized by the three 
parameters n s , A^,(/co), and r. 

In the following, we first consider constraints on n s as- 
suming r = 0, and then on both n s and r. We then 
compare the constraints in the n s -r plane to predictions 
of inflationary models. Constraints on the scale depen- 
dence of the spectral index (dn s /dhik) are considered in 
the companion paper H12. 

6.9.1. Constraints on the Scalar Spectral Index 

Inflation is a nearly time-translation invariant state, 
however this invariance must be broken for inflation to 
eventually come to an end. The wavelength of perturba- 
tions depends solely on the time that they were produced, 
thus a time-translation invariant universe would produce 
scale-invariant perturbations (n s = l). 34 The predic- 
tion that inflation should be nearly, but not fully, time- 

34 Scale invariance here means that the contribution to the rms 
density fluctuation from a logarithmic interval in k, at the time 
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Fig. 8. — The data strongly prefer departures from a scale invariant primordial power spectrum (n s < 1), as predicted by inflation. 
Left panel: The marginalized one-dimensional constraints on n s for the standard ACDM model (with r = 0) using several datasets. SPT 
data tightens the constraint on n s relative to WMAP7 alone. Adding BAO data further tightens this constraint and leads to a preference 
for lower values of n s , while adding Ho has little effect. Right panel: The one-dimensional marginalized constraints on n s from the 
SPT+ WMAP7+Hq+BAO dataset given three different models. Plotted are ACDM (black solid line), ACDM+Sm„ (purple dashed 
line) as a typical case for extensions affecting the late-time universe, and ACDM+AT e ff (green dot-dashed line) as a typical case for 
extensions affecting the Silk damping scale. Of the extensions considered here, only those that affect the damping tail - in this case by 
varying neutrino species - causes noticeable movement towards n s = 1. We note that in all cases the data robustly prefer a scale-dependent 
spectrum with n a < 1. 



translation invariant gives rise to the prediction that n s 
should deviate slightly from unity ( Dodelson et al.|1997 ). 

Because of the special status of n s = 1 , and because we 
generally expect a departure from n s = 1 for inflation- 
ary models, detecting this departure is of great interest. 
Kll combined data from SPT and WMAP7 to measure 
a 3.0 a preference for n s < 1 in a ACDM model, with 
n s = 0.966 ± 0.011. We show our constraints on n s 
for the ACDM model from several combined datasets in 
the left panel of Figure [8j All datasets strongly prefer 
n s < 1. 

Using SPT+ WMAP1 data, we find 



= 0.9623 ±0.0097. 



(26) 



For this dataset, we find P(n s > 1) = 4 x 10~ 5 , a 3.9 cr 
departure from n s = 1; n s < 1 is heavily favored. 

Including BAO data substantially shifts and tightens 
the constraints on n s , as can be seen in Figure [8] The 
BAO distance measure r s /Dy depends on Q\, breaking 
the partial degeneracy between Q\ and n s in the CMB 
power spectrum. The BAO preference for lower Q\ pulls 
the central value of n s down to n s — 0.9515±0.0082. Us- 
ing a high-temperature MCMC, we measure the proba- 
bility for n s to exceed one to be 1.1 x 10 -9 , corresponding 
to a 6.0cr detection of n s < 1. 

Including Hq data has a smaller effect on the n s con- 
straint than BAO, slightly disfavoring low-n s values as 
seen in Figure [HJ The mechanism for the improvement is 
the same as for BAO, however, the CMB and H datasets 
individually prefer similar values of Q\. Thus the two 
datasets tighten the n s constraint around the CMB-only 
value, leading to n s — 0.9638 ± 0.0090. Using the com- 

when k = aH, is independent of k. Here a(t) is the scale factor 
and H = a/a is the Hubble parameter. 



bined CMB+f/o dataset, we measure the probability for 
n s to exceed one to be 3.1 x 10~ 5 , corresponding to a 
4.0(7 preference for n s < 1. 

As expected, combining CMB with both BAO and H 
data nudges the constraint on n s up slightly from the 
CMB+BAO constraint to n s = 0.9538 ± 0.0081, thus 
weakening the preference for n s < 1 slightly from 6.0 to 
5.7cr. 

In summary, regardless of which datasets we use, the 
data strongly prefer n s < 1 in the ACDM model. 

The importance of detecting a departure from scale 
invariance leads us to review our modeling assumptions. 
Specifically, are there extensions to the standard ACDM 
model that could reconcile the data with a scale-invariant 
spectrum, n s = 1? We answer this question by cal- 
culating the n s constraints from the CMB+if +BAO 
dataset for several physically motivated ACDM model 
extensions. 

We consider two classes of model extensions: those 
that can affect the slope of the CMB damping tail, and 
those that cannot. As a representative case of the first 
class of extensions, we consider ACDM+A^ff, in which 
the number of relativistic species is allowed to vary. As 
an example of the second class of extensions, we consider 
massive neutrinos ACDM+Sm^ (with N e g fixed at its 
fiducial value of 3.046). These example extensions as 
well as several others are explored in considerable detail 
in H12. 

Of the extensions considered, only models that can 
affect the slope of the damping tail significantly in- 
crease the likelihood of n s — 1. The results of this 
test are displayed in the right panel of Figure |8j where 
we show the marginalized constraints on n s from the 
CMB+BAO+iJo dataset. Even the N e g extension does a 
poor job reconciling the data with a scale-invariant spec- 
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Fig. 9. — This figure highlights the contributions of the SPT 
data to constraints on the tensor-to-scalar ratio, r. We show 
four datasets: WMAP7 (red dashed line), WMAP7+BAO 
(orange dot-dashed line), SPT+ WMAP7 (black solid line), 
and SPT+ WMAP7+BAO (black dotted line). Note that the 
WMAP7+BAO and SPT+ WMAP7 likelihood functions are nearly 
identical. SPT data tightens the r constraint significantly, regard- 
less of whether BAO data are included. While adding low-redshift 
Ho measurements has minimal effect on the constraints on r (not 
shown), adding low-redshift information from BAO tightens con- 
straints on r considerably. SPT+WMAP7 constrains r < 0.18 
(95% C.L.), while adding low-redshift BAO measurements tight- 
ens the constraint to r < 0.11 (95% C.L.). 

trum; the cumulative probability for n s > 1 is 6.1 x 10~ 3 , 
constraining n s to be less than one at 2.5 a. The pre- 
ferred value for N e g is far greater than the nominal value 
of 3.046 in the limited allowed parameter space. 

We conclude that the data robustly prefer a scale- 
dependent spectrum with n s < 1. 

6.9.2. Constraints on Tensor Perturbations 

The last inflationary parameter we consider is the 
tensor-to-scalar ratio r. Because r is related to the en- 
ergy scale of inflation, 35 a detection of r would provide 
an extremely interesting window onto the early universe. 
We first consider the marginalized constraints on r for 
the ACDM+r model, shown in Table [4] and Figure [9j 
then move on to a comparison with inflationary models 
in the n s -r plane in § |6.9.3| and Figure [10] 

One can think of the r measurement in the follow- 
ing way. The CMB power spectrum is first measured 
at £ > 60 (where tensor perturbations are negligible) to 
determine the ACDM model parameters, thus determin- 
ing the scalar contributions to the power spectrum. This 
scalar-only spectrum is then extrapolated to low £; any 
excess power observed is due to tensor perturbations. In 
this way, although the SPT data presented here do not 
directly measure power that could be from gravitational 
waves, by pinning down other model parameters, the ex- 
trapolation of the scalar power to large scales is more 

35 The tensor-to-scalar ratio r is proportional to the inflaton 
potential V(rfi) and the energy scale of inflati on is proportional to 
V(<j>) 1 / i I Kinney|2003 Baumann et al.|2009 l. 



precise. 

Even with the parameters that determine the scalar 
power spectrum perfectly known, there would still be 
significan t uncertainty in the v alue of r due to cosmic 



variance. |Knox fc Turner ( 1994 ) showed that for different 
realizations of a universe where r = 0, 50% of cosmic 
variance limited full sky temperature surveys will be able 
to place a limit of r < 0.1 at 95% confidence. As we will 
see, the results we present here approach that limit. 
WMAP7 data alone have been used to constrain r < 



0.36 at 95% confidence in the ACDM+7- model (Larson 
et al.|20lT ) . Prior to the current SPT analysis, t he tight 



est published constraint on r was reported by | Sanchez] 
|et al. (2012), who used the combination of datasets from 
TJOm^CMASS, WMAP7, and Kll to constrain r < 0.16 
at 95% confidence. 

The bandpowers presented here lead to a significant re- 
duction in the upper limit on r. These measurements are 
summarized in Table|4j As shown in Figure[9j SPT band- 
powers tighten the constraint on r regardless of whether 
low-redshift information from BAO or Hq is included. 
Using the CMB datasets, we measure 



r < 0.18 (95% C.L.) . 



(27) 



This limit remains unchanged for the CMB-|-i?o datasets, 
while from the CMB+BAO datasets we measure r < 0.11 
(95% C.L.). As with constraints on n s , the improvement 
with the addition of BAO comes from the improved BAO 
constraints on fljy breaking a partial three-way degener- 
acy between ft\, n s , and r. In contrast, adding in the H 
measurement to CMB data does not significantly tighten 
the CMB-only result because the Ho measurement sharp- 
ens up the Q\ distribution around the values that allow 
for larger r. It is worth noting that the SPT bandpowers 
tighten the constraint on r even in the presence of BAO 
data; this can be seen in the tighter r constraint from the 
WMAP7+BAO to the SPT+ WMAP7+BAO datasets in 
Figure [9j 

Using the combination of the CMB+BAO+-ff 
datasets, we measure 



r < 0.11 (95% C.L.) ; 



(28) 



this limit is unchanged from the CMB+BAO constraint. 
The Hq and BAO data pull in opposite directions on the 
CMB-only data in the r s /Dy - Hq plane, however the 
statistical weight of the BAO measurements dominate 
the constraint in the combined dataset. 

With the addition of low-redshift information from 
the BAO measurement, the constraints on tensor per- 
turbations approach the theoretical limit of what can 
be achieved with temperature anisotropy alone. With 
the combination of these data, the cosmological param- 
eters controlling the scalar perturbations are now suf- 
ficiently well known that they do not significantly de- 
grade the limit on r. This is evident in Figure [TO} where 
adding BAO removes most of the degeneracy between 
n, and r. Further improvements now await precision 
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C.L.) from the BICEP experiment ( |Chiang et aT|2010| ). 
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Fig. 10. — We compare the constraints on the ACDM+r model with predictions from models of inflation in the n a — r plane. We show 
the two-dimensional constraints on r and n s as colored contours at the 68% and 95% confidence levels for three datasets: WMAP7 (grey 
contours), CMB (red contours), and CMB+iio+BAO (blue contours). Adding the SPT bandpowers partially breaks the degeneracy 
between n s and r in the WMAP7 constraint, which can be seen clearly moving between the grey and red contours. Plotted over the 
constraint contours are predictions for several models of inflati on. We restrict our com parison with model predictions to the simplest cases 
of slow-roll inflation due to a single scalar field as reviewed in |Baumann ct al. (2009}. 

Solid black line: The predictions of exponential inflation (V((f>) oc exp ^yj 16ir<f> 2 / (p TVfpj)J ) lie on this line. In exponential inflation, 
increasing p moves the prediction towards the Harrison-Zel'dovich-Pccblcs point n s = 1, r = 0. 

Black lines with colored circles: The predictions of chaotic inflation models (V(4>) oc ((f>/fj,) p , p > 0) for five different values of p lie on 
the corresponding line. The predictions in the r - n s plane are a function of TV, where TV is expected to be in the range TV £ [50, 60]. 
Purple region: The predictions of large-field hill-top inflation models (V(<j>) oc 1 — (0//x) 2 ) lie within the colored region, which is shown 
for TV £ [50, 60]. 



We now turn to a comparison of model predictions 
with data constraints in th e n s -r plane. This compar- 
ison is illustrated in Figure 



10 where we show the two- 



dimensional marginalized constraints from three combi- 
nations of data with predictions from simple models of 
inflation over-plotted. First, we note that the confidence 
contours for the CMB-only case in Fig. [10] show the ex- 
pected positive correlation between n s and r. Essen- 
tially, the suppression of large-scale power when increas- 
ing n s can be countered by adding extra large-scale power 
sourced by tensors. The SPT data disfavor large values 
of n s (and hence r) , significantly reducing the degeneracy 
between these two parameters. With BAO data added to 
SPT+ WMAP7, the n s — r correlation nearly disappears. 
Adding Hq data has little effect on constraints from the 
CMB or CMB+BAO datasets, removing only the small- 
est allowed values of n s in both cases. As mentioned 
above, we are approaching the cosmic variance limit for 
the temperature anisotropy on measuring r - at which 



point improved knowledge of the six ACDM parameters 
no longer translates into better limits on r. 

We restrict the model comparisons to the simplest 
cases of single-fi e ld, sl ow-roll inflation, as reviewed in 



Baumann 
acterized 



ct al. 



(|2009|). Models can be broadly char- 
how much the inflaton field 6 



according to 

changes from the time perturbations on observably large 
scales were being produced until the end of inflation; this 
change in 4> defined at Acf>. Models in which Acf> is larger 
than the Planck mass (Mp\) are classified as "large- field" 
models, while those in which A<p < Mp\ are classified as 
"small-field" models. The dividing line between the two 
cases corresponds to r = 0.01. 

Here we look at large-field inflation models, consid- 
ering several forms of the inflaton effective potential: 
large-field chaotic inflation models (V((p) oc (0//x) p , 
p > 0), large-field hill-top inflation models (V((f>) oc 
1 — (<f)/[i) 2 ), and exponential inflation models [V(<f>) oc 

exp \y/mw(j> 2 /( P Mgj]). 
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Large-field chaotic models have the fewest free parame- 
ters, and we discuss them first. Given a choice of p, these 
models have just one free parameter, and this parameter 
is highly constrained by the requirement of reproduc- 
ing the well-known amplitude of the scalar perturbation 
spectrum. Thus, these models make fairly localized pre- 
dictions in the n s — r plane. The uncertainty in these 
predictions is dominated by the details of the end of infla- 
tion which are not specified by V((f>) but, instead, depend 
on the coupling of the inflaton field <j> to other fields. This 
uncertainty can be captured by the parameter N where 
e N gives the increase in the scale factor between the time 
when the observable scale leaves the horizon and the end 
of inflation. 36 Assuming a standard slow-roll inflation 
scenario, 37 A^ is expec ted to lie in the range 50 to 60 



(Liddle & Leach 2003). The spread in values is dom- 
inated by uncertainty in how much the energy density 
drops between the end of inflation and reheating, though 
this range can be extended in either direction by modifi- 
cations to the standard thermal history. In Figure 10 we 
consider several large-field chaotic models, each with a 
different value of p, and indicate the predictions of each 
model as A" varies between 50 and 60. 

The A(j) 4 (p = 4) m odel was ruled out by earlier data 
( Komatsu et al.|2009 |, given the expected range of N. It 
possibly could have been saved with a non-standard post- 
inflation thermal history (designed to make N very large, 
thus m ovi ng the prediction towards r = and n s =l in 
Figure 10), but such a maneuver no longer works; the 
CMB+BAO dataset excludes the X(f> 4 model with greater 
than 95% confidence for all values of N. Similarly, the 
<f) 3 (p — 3) model is excluded with greater than 95% 
confidence by the CMB or CMB+ifg datasets given the 
expected range of N, and excluded with greater than 
95% confidence by the CMB+BAO datasets regardless 
of N. 

While the m 2 4> 2 (p = 2) model was consistent with 



previ ous constraints ( Keisler et ah] 2011 Sanchez et al. 
2012 1, the current combinations ot CMB with BAD data 



place a tight upper limit on r and disfavor the m 2 </> 2 
model, which produces predictions that fall at the edge 
of the 95% confidence contour. This model is allowed by 
CMB or CMB+TJo- Models with smaller values of p are 



consistent with the data, as shown in Figure 10 

The exponential models lead to n s and r predictions 
that are independent of scale and therefore independent 
of A". The predictions, as p varies, form a line in the 
n s — r plane. This whole class of models is allowed at 
95% confidence for a range of p by the CMB+H data, 
and is excluded (> 95% C.L.) by the CMB+BAO and 
CMB+ffo+BAO data. The potential in hill-top mod- 
els has the shape of symmetry-breaking potentials that 
drive <f> away from the origin. The generic form of the 
hill-top potential is V(cj)) oc 1 — ((f>/p) p . Such models, for 
a fixed p, have two free parameters, the proportionality 
constant and p. The proportionality constant is set by 

36 A note of clarification: N says nothing about the total increase 
in the scale factor between the beginning and end of inflation, which 
is expected to be much larger. 

37 The assumption here is that inflation stops by the end of slow 
roll and is followed by the field oscillating in an approximately 
quadratic potential near the minimum. The universe eventually 
reheats to a density greater than that during BBN, and then the 
standard thermal history ensues. 



the amplitude of scalar fluctuations A|j, while /i is con- 
strained by n s and r. For p < 2, hill-top models have the 
behavior of large-field models for the range of n s allowed 
by the data; the predictio n of the p = 2 case is shown as 
a purple region in Figure |10[ and is consistent with the 
data. 

Finally, we consider small-field inflation models. Ex- 
amples of small-field potentials include hill-top poten- 
tials with p > 2 (fo r some values of p), the Co leman- 
Weinberg potential (Coleman & Weinberg 1973) (with 
suitabl y adjusted paramete rs), and warped D-brane in- 
flation ( Kachru et al.|[2003 ). Because small-field models 
predict r < 0.01, the SPT data have no constraining 
power on these models through limits on r. All of these 
models are consistent with the data, so long as they are 
consistent with the limits on n s . 

Constraints on the scale dependence of the spectral 
index (dn s /dh\k) are considered in the companion paper 
H12. 

7. CONCLUSION 

We present a measurement of the CMB temperature 
power spectrum from 2540 deg 2 of sky observed with the 
South Pole Telescope (SPT). These are the first CMB 
power spectrum results reported for the full SPT-SZ sur- 
vey, which encompasses three times the ar ea used in 
previous SPT power spectrum analyses (Kll, Reichardt 
et al. (2012b)) The bandpowers cover the third to ninth 
acoustic peaks (650 < £ < 3000) with sample-variance- 
limitcd precision at £ < 2900. This measurement repre- 
sents a significant advance over previous measurement s 
of the damping tail by ACBAR (Reichardt et al. ||2009j, 
QUaD ([Brown et . al.|[2009| [Friedman et a l. 2009), ACT 



( |L>as et al ||2Ullri and SPT ( |Keisler et al.p01l| ). 

We find the SPT bandpowers are well fit by a spa- 
tially flat ACDM cosmology with gravitational lensing by 
large-scale structure. We use this SPT measurement to 
extend the dynamic range probed by the WMAP power 
spectrum, thus tightening parameter constraints in the 
six-parameter ACDM model. With the exception of the 
optical depth r which is constrained by the large-scale 
polarization data from WMAP7, adding the full survey 
SPT bandpowers significantly improves measurements of 
all ACDM parameters. Most notably, the measurement 
of the angular sound horizon, 8 S , tightens by a factor of 
2.7 due to the number of acoustic peaks detected at high 
signal-to-noise. Uncertainties on the other four parame- 
ters are reduced by a factor of ~ 1.4. The combination 
of SPT and WMAP7 data is used to constrain n s < 1 at 
3.9 a. 

We examine constraints on three extensions to the 
ACDM model. We first allow for a rescaling of the 
gravitational lensing potential by a parameter A^. Us- 
ing CMB data, we rule out the no-lensing hypothesis 
(Al < 0) at 8.1a, the most significant detection to date 
using the CMB alone, and measure a lensing amplitude, 



.4 



+0.15(+0.30) 



l = 0.86I^3|Io 25j (68% and 95% C.L.), consistent 

with the ACDM expectations. We expect the lensing de- 
tection significance to triple in a future analysis of the 
full SPT-SZ survey using an optimized four-point lens- 
ing estimator, similar to the one appli ed to one fifth of 
the survey by van Engelen et al. (2012). 

Second, we allow non-zero mean curvature of the ob- 
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servable universe. The low-redshift information encoded 
in the CMB by gravitational lensing helps to improve 
constraints on the mean curvature of the observable uni- 
verse for the ACDM+ri fc model. Using the CMB alone, 
we measure fi^ = — 0.003+q oil) which is consistent with 
a flat universe. Models without dark energy are ruled 
out at 5.4 a. 

Finally, we look at constraints on the amplitude of ten- 
sor perturbations. The combination of SPT+WMAP7 is 
used to constrain the tensor-to-scalar ratio to be r < 0.18 
with 95% confidence. 

Adding low-redshift probes of Ho and BAO further 
tightens these constraints. Combining the CMB and Ho 
datasets mildly tightens our parameter constraints, and 
is fully consistent with the CMB constraints. Combining 
the CMB and BAO datasets, on the other hand, leads to 
significant improvements in our parameter constraints. 
Combining all three datasets produces constraints that 
lie close to the CMB+BAO constraints; the combination 
of CMB+i?o+BAO is used to constrain n s < 1 at 5.7 a 
in the ACDM model, measure Cl k = - 0.0059 ± 0.0040 
in the ACDM+fi fe model, and constrain r < 0.11 at 95% 
confidence in the ACDM+r model. This constraint on r 
approaches the theoretical limit of how well tensor per- 
turbations can be constrained from CMB temperature 
anisotropy, r < 0.1 (95% C.L.). We compare these con- 
straints on n s and r to the predictions of single-field in- 
flation models and exclude several models with greater 
than 95% confidence. 

There is some tension between the six datasets in- 
cluded in the CMB+BAO+H combination for a ACDM 
cosmology. However, we assume the uncertainties re- 
ported for each of the datasets are correct and combine 
them to produce many of the results presented here. We 
refer the reader to H12 for a more detailed discussion of 
the consistency of the datasets. 

In this paper, we have focused on the amplitude and 
shape of the primordial power spectrum of scalar and ten- 
sor perturbations, as well as the effects of gravitational 
lensing and curvature. Further cosmological implications 
of the bandpowers from the full SPT-SZ survey, includ- 
ing constraints on the neutrino masses, the dark energy 
equation of state, the primordial Helium abundance, and 



the effective number of neutrino species are explored in 
the companion paper H12. 

The first CMB temperature power spectrum from the 
Planck satellite is expected to be released in 2013. We 
expect that, given the small beam of the SPT relative to 
Planck, the data presented here should remain the most 
precise measurement of small angular scale anisotropy 
for 2200 < I < 3000. The cosmological constraints pre- 
sented here will tighten significantly with the with the 
Planck power spectra. Further improvement will come 
with the addition of polarization information from up- 
coming experiments including Planck, ACT-Pol, Polar- 
bear, and SPT-Pol, the new polarization-sensitive cam- 
era on the SPT. 
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